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Among all forms of routine human activity, the one which produces the strongest gravity-gradient noise in 
interferometric gravitational-wave detectors (e.g. LIGO) is the beginning and end of weight transfer from one 
foot to the other during walking. The beginning and end of weight transfer entail sharp changes (timescale r ~ 
20msec) in the horizontal jerk (first time derivative of acceleration) of a person's center of mass. These jerk pairs, 
occuring about twice per second, will produce gravity-gradient noise in LIGO in the frequency band 2.5Hz < 

/ < 1/(2t) ~ 25Hz with the form ^S h (f) ~ 0.6 x 10" 23 Hz- 1/2 (//10Hz)- 6 (^(ri/lOm) -8 ) ^ ■ Here 
the sum is over all the walking people, n is the distance of the i'th person from the nearest interferometer test 
mass, and we estimate this formula to be accurate to within a factor 3. To ensure that this noise is neglible 
in advanced LIGO interferometers, people should be prevented from coming nearer to the test masses than 
r ~ 10m. A r ~ 10m exclusion zone will also reduce to an acceptable level gravity gradient noise from 
the slamming of a door and the striking of a fist against a wall. The dominant gravity-gradient noise from 
automobiles and other vehicles is probably that from decelerating to rest. To keep this below the sensitivity of 
advanced LIGO interferometers will require keeping vehicles at least 30 meters from all test masses. 

PACS numbers: 04.80.Nn 



I. INTRODUCTION AND SUMMARY 

Time-changing Newtonian gravitational forces, acting on 
the test masses of an interferometric gravitational-wave detec- 
tor (e.g. in LIGO), produce noise. This noise is convention- 
ally called "gravity-gradient noise" because the interferome- 
ter measures the differences in the gravitational forces acting 
on the test masses. In a previous paper, Thorne and Hughes 
[|l|] gave a general overview of gravity-gradient noise and ana- 
lyzed in detail the gravity-gradient noise in LIGO due to den- 
sity fluctuations in the earth caused by ambient seismic waves. 
In this paper we focus on gravity-gradient noise due to human 
activity. 

It has long been recognized that gravitational forces 
from moving humans can produce significant noise in high- 
precision gravitational experiments. Roll, Krotkov and Dicke 
[^] took great care to eliminate such forces in their classic 
Eotvos experiment, and Dicke ^,||] has raised the possibil- 
ity that such noise was a serious factor in Baron Roland von 
Eotvos's original versions of that experiment [Q]. In the early 
years of the LIGO Project, Robert Spero and others made 
rough estimates of the magnitude of human (and other animal) 
gravity-gradient noise in LIGO and the distances to which 
humans (and other animals) should be relegated to control 
it. While these estimates were sufficiently accurate for their 
purposes, no analyses until ours seem to have identified the 
form of routine human activity that will dominate the noise 
(human walking), nor the spectrum of this dominant noise, 
oc f- 6 for 2.5Hz < / < 25Hz. 

A first version of our analysis was carried out in summer 
1995 when the buildings that house LIGO's test masses 
were being designed. Our goal was to make sure that the 



building design would keep humans sufficiently far from the 
test masses, during routine LIGO observations, for human 
gravitational noise to be unimportant. The press of other re- 
search delayed until now our finalizing and publishing this 
analysis. 

In our 1995 document |^], we focused on the gravitational 
effects of a walking person's horizontal center-of-mass mo- 
tion, and we based our analysis on two force-plate experi- 
ments taken from the Biokinesiology (motion-in-biological- 
systems) literature; cf. Sec. [I A below. For the present pa- 



per, we have augmented our center-of-mass-motion data base 
with new force-plate measurements on three different persons; 
and using data from the Biokinesiology literature, we have 
verified our original guess that the motion of a walking per- 
son's limbs produces gravitational noi se sm all compared to 
that from center-of-mass motion (Sec. II C). This extended 
analysis has not changed significantly any of our 1995 con- 
clusions. 

Our analysis (Sec. [|) produces the following estimate for 
gravity-gradient noise in LIGO due to walking people: 



0.6 x 10~ 23 /10Hz\ 



VHz" \ f J 
at 2.5Hz < / < 25Hz . 



1/2 



(1) 



Here Sh(f) is the spectral density of the interferometer's out- 
put gravitational-wave signal h = AL/L (with L = 4km 
the interferometer arm length and AL the difference in arm 
lengths, which fluctuates due to the gravitational forces from 
the walking people); also r"j is the distance from person i to 
the nearest interferometer test mass. 
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We believe this estimate to be accurate to within a factor 3. 
A factor ~ 2 uncertainty arises from the angular location and 
direction of motion of each person [the factor a in Eq. ( |l2| ) 
below], and another factor that on occasion may be as large as 
~ 2 arises from the gravitational forces of compressional and 
shear waves in the floor and ground, produced by the person's 
walking (Sec. [IB). Somewhat smaller than this are uncer- 
tainties due to variations in the gait (walking) pattern from 
one person to anoth er, and for each person, from one step to 
another (Sec. [I A 2 ). Adding our two factor ~ 2 uncertainties 
in quadrature, we get our net factor ~ 3 uncertainty. 

The 1/rf dependence of the noise ([!]) results from the fact 
that it is produced by changes in the distance to the person's 
center of mass, and thus by the person's changing dipole grav- 
itational field. The 1/ f 6 dependence results from two facts: 
(i) The gravitational force produces a test-mass acceleration, 
which means a second time derivative of AL and thence a 
second time derivative of h, and thence a l// 2 in the am- 
plitude spectrum \fSh- (ii) The fourth time derivative of the 
center-of-mass position is the lowest-order derivative that has 
a delta-function-like behavior on timescales short enough to 
produce noise at frequences / ~ 10 Hz; and those four time 
derivatives produce an additional l// 4 in the spectrum. 

Figure [l] shows the noise ([!]) for a single person at vari- 
ous distances from the nearest test mass. For comparison 
we also show (i) the benchmark noise curve for a broad-band 
"advanced" interferometer (which might operate in LIGO in 
the ~ 2010 time frame) (ii) the standard quantum limit 
(SQL) for a LIGO interferometer with 1 tonne test masses, 
and (iii) the estimated seismic gravity-gradient noise from am- 
bient earth motions [[j]] at "quiet times" (assuming the "stan- 
dard LIGO seismic spectrum") and at "very quiet times" (as- 
suming a seismic spectrum 10 times lower than the standard 
one - a level that might occur during wind-free nights). 

Figure [l] shows that a single person walking at a distance 
of 5m from a test mass could significantly increase the noise 
in an advanced interferometer, and several people would be 
correspondingly more serious. 

The LIGO corner building (the only one with extensive hu- 
man activity) has been designed to keep people at least 10 
meters from all test masses during normal operations. This 
provides an adequate safety factor for advanced interferome- 
ters; if the noise is three times as large as our estimate, then 
10 people at 10 meters distance would increase an advanced 
interferometer's noise by only a few tens of percent near the 
most sensitive frequency, 10Hz. 

It is conceivable that in ~ 2010 or later interferometers will 
be operated in LIGO with good performance at frequencies 
below 10 Hz — e.g., as low as ~ 3 Hz. (Adalberto Giazzotto 
and colleagues of the VIRGO Project have developed seis- 
mic isolation systems that can go down to frequencies as low 
as 3 Hz [gp.) Figure [j] indicates that, if such interferometers 
are ever operated in LIGO, it will be necessary to expand the 
people-free zone around each corner test mass. 

The LIGO end and mid-station buildings (which have lit- 
tle human activity) are designed to keep all humans at least 
5 meters from the end test masses during normal operations. 
This provides an adequate safety margin for the first and en- 
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FIG. 1. The predicted spectrum of noise in LIGO's 4km-long in- 
terferometers due to human gravity gradients (dark strips) and seis- 
mic gravity gradients (light strips), compared with the benchmark 
noise curve for advanced LIGO interferometers and with the standard 
quantum limit (SQL) in interferometers with 1 tonne test masses. 
The thickness of the strips indicates the estimated uncertainties in 
our analysis. 



hanced LIGO interferometers, which will operate in the early 
and mid 2000's; but when advanced interferometers begin to 
operate, it will be necessary to expand the people-free zone, 
most especially at the end of each end building. 

Robert Spero's early estimates of human gravity gradient 
noise focused on the slamming of a door or the striking of a 
fist against a wall [g]. Since this is more readily suppressed 
(by warnings and viscous door stops) than human walking, we 
have regarded it as a less serious and pervasive noise source. 
However, whenever a door slams or a fist strikes, the mag- 
nitude of the resulting gravitational "signal" in an advanced 
LIGO interferometer will be comparable to that from people 
walking: 

In Sec. |l| by a variant of Spero's analysis we derive the 
following expression for the Fourier transform of the signal 
h[t) produced by a mass M striking a building wall and com- 
ing suddenly to rest: 



\h\ = 



GM\aAv\ 
Lr 3 (27r/) 4 



(2) 



Here Av is the object's sudden change of speed, r is its dis- 
tance from the nearest interferometer test mass, L = 4km is 
the interferometer arm length, and a is a coefficient in the 
range —2 < a < +2 that depends on object's angular loca- 
tion. Following Spero, we show that, with optimal signal pro- 
cessing, this "signal" would produce the following amplitude 
signal-to-noise ratio in an advanced LIGO interferometer: 



10m 



S_ ^ ( MAv 
N \ 5kg m/s y \ r 



(3) 



Here our fiducial value, 5 kg m/s, for MAv corresponds to a 
5 kg door slamming shut at a speed of 1 m/s, or a 2 kg fist and 
forearm striking a wall at a speed of 2.5 m/s. 
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Thus, slamming doors and striking fists, like walking peo- 
ple, must be kept at distances r > 10m from the test masses 
of advanced LIGO interferometers. 

A third type of human activity that can produce strong grav- 
itational noise is the motion of automobiles and other vehicles. 
In Sec. |lli| we argue that the dominant vehicle gravity-gradient 
noise, in the critical frequency region / ~ 10 Hz, is due to a 
vehicle's sudden change of acceleration when it comes to rest, 
e.g. when parking. If Aa is the vehicle's change of acceler- 
ation, M is its mass, and r is its distance from the nearest 
interferometer test mass, then in an advanced LIGO interfer- 
ometer the vehicle will produce a "signal" h(t) with Fourier 
transform in our frequency band 



\h\ = 



GM\aAa\ 
Lr 3 (2nf) 5 ' 



(4) 



Here a is the same angle-dependent coefficient with range 



-2 < 



< 



-2 as appears in the door/fist analysis. With op- 



timal signal processing, this "signal" would produce the fol- 
lowing amplitude signal-to-noise ratio in an advanced LIGO 
interferometer: 



--if— 1 
N V 2tonne / 



O.65 / \ r 



^30m\ 3 



(5) 



Here g — 9.8ms~ 2 is the acceleration of gravity and we have 
set I a I to a representative value, \/2. 

The LIGO service roads come no closer than 40 meters to a 
LIGO corner test mass, but they approach to within 15 meters 
of the corner- and mid-station test masses. 

The product M|A<x| = 2tonne x 0.6g used in Eq. (Q) is in 
the upper range of what one might expect for stopping vehi- 
cles. Two tonnes is a modest vehicle mass; O.65 is the decel- 
eration at which a vehicle begins to skid on dry asphalt. Thus, 
when advanced interferometers are operating in LIGO it will 
be necessary to increase the radius of the vehicle-free zone at 
the corner and mid stations to > 30m. 

The body of this paper is organized as follows: At the be- 
ginning of Sec. [| we briefly explain why human walking is 
the dominant source of human gravity gradient noise in inter- 
ferometers, and why, at the frequencies of interest, the noise 
comes pred ominantly from sudden changes in motion. Then 
in Sec. |5 A| we compute the gravity-gradient noise in an inter- 
ferometer due to a person's center-of-mass motion, and sum 
over a population of people to get Eq. ([j]) (discussed above). 
Momentum conservation implies that any sudden change in 
a person's center-of-mass motion will produce a correspond- 
ing sudden displa cem ent of the floor and the ground beneath 
the floor; in Sec. [IB and an Appendix we show that gravity 
gradient noise from this floor/ground motion will not cancel 
that from the person, but can, on occasion, cancel as much as 
half of it (thereby intro ducing a factor ~ 2 uncertainty in the 
net noise). In Sec. |l C| we compute the gravity-gradient noise 
due to the motions of a person's limbs (including, most im- 
portantly, the sudden changes of motion when a heel strikes 
the floor). We show that the limb motion produces noise that 
is smaller by a factor ~ 0.1(10m/r)(//10Hz) than the noise 



from center-of-mass motion; here r is the distance between 
the person and the nearest test mass. In Sec. |l|we analyze the 
gravity-gradient noise due to the sudden stopping of a moving 
mass — a slamming door, a fist striking a wall, or a parking 
vehicle — arriving at Eqs. (|j) and (Q), discussed above. In Sec. 
[V we make some concluding remarks. 



II. HUMAN WALKING 



Consider a person (or vehicle) moving in the vicinity of a 
LIGO test mass. Denote by $ the person's Newtonian grav- 
itational potential in the test-mass vicinity, and by x'{t) and 



r'(t) 



the vector and distance from the person's cen- 



ter of mass ("CofM") to the test mass at time t. Expand the 
Newtonian potential in multipole moments around the per- 
son's (moving) CofM: 



$ = 



GM 3 GX jk 



x j x k 



(6) 



Here Tjk(t), the person's quadrupole moment relative to his 
CofM, is the means by which his moving limbs produce 
gravity-gradient noise, repeated indices (j and k) are to be 
summed, and we use Cartesian coordinates so it doesn't mat- 
ter whether tensor indices are up or down. 

In Sec. p A| we examine the first (monopolar) term in Eq. 
(^). This is the dom inant gravitational effect of the CofM 
motion. In Sec. QC we examine the second (quadrupolar) 
term — including its time dependence due both to limb mo- 
tions and to center-of-mass motions — and we show that its 
influence on an interferometer's noise is small compared to 
that of the monopolar term. 

Before presenting these analyses, it may be useful to com- 
ment on the importance of sudden changes (jerkiness) in the 
human (or vehicular) motion. 

Smooth (non-jerky) motion produces gravitational forces 
that are concentrated at frequencies / ~ 1 Hz, well below 
those, / > 10 Hz, of concern for the interferometers. For ex- 
ample, the period of the normal human gait cycle (two steps, 
one left and one right) is about 1 sec (frequency 1 Hz); and 
an autumobile moving at speed 30km/hr at a distance 15m 
from a test mass travels through an angle 90° as seen by the 
test mass in about 2 sec (frequency ~ 0.5 Hz). If the motion 
is sufficiently smooth, the Fourier transform of such motions 
will fall off with frequency exponentially, becoming totally 
negligible at / ~ 10 Hz. 

By contrast, if the rt'th time derivative of the motion 
changes significantly on a timescale r < 50msec, then the 
n+ l'th time derivative will have a sharp, delta-function-like 
peak with time width r, and that will produce a Fourier trans- 
form of the motion that falls off as 1/ f n+1 up to frequencies 
/ ~ 0.5/r > 10Hz. Such a power-law fall-off can produce 
much larger gravity-gradient noise near 10 Hz than the expo- 
nential fall-off caused by smooth motion. 
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A. Motion of Center of Mass 

In this section we examine the gravitational noise produced 
by the monopolar (center-of-mass) term in Eq. (H). 



1. General formulas for noise 

We focus on a time duration of one gait cycle ~ lsec, cen- 
tered on time t = 0. We introduce Cartesian coordinates Xj 
attached to the floor, with origin at the location of the person's 
CofM at time t = 0; and we denote by £ (t) the motion of the 
CofM relative to this origin (so £ = at t = 0). Then the 
vector and distance x' and r' from the person's CofM to the 
test mass are 



-el 



2r 



(7) 



Here x and r are the vector and distance from the origin of 
coordinates to the test mass, and h = x/r is the unit vector 
pointing from the origin to the test mass. Correspondingly, we 
can write the center-of-mass piece of the gravitational poten- 
tial (||) as 



GM 



CofM 



GM GMZjUj 



(8) 



3 GMfatk - l\j\ 2 S jk ) nj n k 
2 r 3 

The first term (monopolar about our fixed center of coor- 
dinates) is constant in time and thus can produce no gravity- 
gradient noise, so we shall ignore it. The second term (dipolar 
about our fixed center of coordinates) produces the gravity- 
gradient noise via sudden changes of the CofM location £. 
The third term (quadrupolar about our fixed center of coordi- 
nates) is smaller than the second, dipolar term by ~ ^ 
0.7m/10m ~ 0.1, and its gravity-gradient noise is corre- 
spondingly smaller in our frequency band, so we shall ignore 
it. The gravitational acceleration of the test mass, produced 
by the person's CofM motion, is then minus the gradient of 
the second, dipolar term: 

GM 

9 j = — g-(£j - injUk^k) . (9) 

The interferometer has four test masses labeled A = 
1,2,3,4, each of which experiences an acceleration of the 
form ([})) due to the CofM motion. The resulting output sig- 
nal of the interferometer, h(t) = AL(t)/L, has a second time 
derivative given by the sum of these accelerations projected 
onto unit vectors rh A that point along the interferometer arms 
as shown in Fig. |[ 



d 2 h _ 1 ^ 



JA 

A 

GM 



m A 



fn A - 3{h A ■ m A )(h A ■ fj 



(10) 
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FIG. 2. Geometry of the LIGO test masses (solid squares) and the 
location (large dot) of the center of mass of a person at time t — 0. 



Since the interferometer's laser beams are horizontal, the 
vectors fh A are all horizontal; and since the person's CofM 
is at the same height as the test masses to within < 1 meter 
and the person is at a distance ~ 10 m from the nearest test 
mass, the vertical component of n A is < 0.1 of the horizontal 
component. This means [cf. Eq. (|Tc|)] that the vertical compo- 
nent of £ produces significantly less gravitational noise than 
the horizontal component: its contribution to \fSh is less by a 
factor < 0.1|£ v |/|£h| ~ 0.3. (Here the tilde denotes a Fourier 
transform, v and h denote vertical and horizontal components, 
and we have used | £ v | / 1 £h | ~ 3 at J ~ 1 Hz as inferred from 
force-plate measurements; Table | below.) On this basis, we 
shall ignore the vertical component of the CofM motion and 
approximate £ and n A as purely horizontal. 

It will be convenient below to rewrite Eq. ( |To| ) in the form 



GMj 



(11) 



where r is the distance from the center of coordinates to the 
nearest test mass, £(t) (scalar, not vector) is the distance the 
CofM has traveled since t = (£ = — 1£| for negative times 
and +|£| for positive), and a is a dimensionless coefficient 
given by 



E - 

A W 



£ • m A - 3(rU • rh A )(h A ■ £) 



(12) 



Here £ = £/£ is the unit vector along the CofM's direction of 
motion, which we regard as constant during one gait cycle. 

For people near an end mass (mass A = 3 or 4 in Fig. ||), 
only that mass gives a significant contribution to the coeffi- 
cient a, and it is easy to verify that a ranges from -2 to +2, 
depending on the angular location of the person. The extremal 
values ±2 are reached when the person's CofM is along the 
interferometer arm and the person is moving toward or away 
from the test mass (£, fh A , and n A all parallel or antiparallel). 

For people in LIGO's corner station, both corner masses 
contribute strongly to a. This is because the distance between 
the two corner masses, I ~ 5m, is comparable to the distance 
between the person and the nearest test mass, r ~ 10m. A 
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FIG. 3. The forward, horizontal force exerted on the floor by a 
woman weighing 73 kg, as measured using dual (two-feet) force 
plates by Earnest L. Bontrager [ jl2| ] (data set 3506a5). Upper panel: 
The force as a function of time exerted by each foot. Lower panel: 
The sum of the forces from the two feet; the full gait cycle is divided 
into two half cycles, A and B. 



straightforward numerical exploration shows, in this case, that 
a ranges from about -2.2 to about +2.2 depending on the per- 
son's location. 

A representative value for \a\, which we use in our final 
noise estimates [e.g., Eq. (0) above], is | a [representative = V%- 



Eq. ( p"lj ) implies that the noise spectrum yj ~ShJj) will scale 
with frequency as \£(f)\/ f 2 , where £ is the Fourier transform 
of the distance traveled £(t); and the discussion just before 
the beginning of this section implies that £(/) will be gov- 
erned primarily by the lowest order derivative of £(t) that has 
sudden changes on a timescale r < 50 msec. To identify that 
derivative and the details of the sudden changes, we rely on 
experimental data from the field of Biokinesiology (the study 
of motion in biological systems). 

Since £ (t) is distance traveled by the CofM, the only way 
that it or its derivatives can change jerkily is by the appli- 
cation of a sharply changing, horizontal external force. The 
only such force, as a person walks, is the horizontal force 
of the floor on the person's feet. The negative of that force 
(the horizontal force F(t) of the feet on the floor) is measured 
by biokinesiologists, using force plates (pp. 414^418 of Ref. 
[]1C}]; Sec. 4.2 of Ref. [[Til]). By momentum conservation, this 
measured force is F — —McP^/dt 2 , and correspondingly the 
gravitational noise is related to the measured horizontal force 
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2. Force-plate measurements 

Figure |] shows the measured horizontal force F(t) exerted 
on the floor by a woman weighing 73 kg during a full gait 
cycle. These data were obtained as follows: 

Our colleague, Earnest L. Bontrager placed two force plates 
in the floor of his laboratory at Rancho Los Amigos Medical 
Center, Downey, California. The force plates were so located 
that in normal walking a person will encounter them during 
one gait cycle, with the right foot landing on the first plate 
and then the left foot on the second. Each plate was equipped 
with piezzo-electric transducers that measured the foot's for- 
ward horizontal ("progressive") force, its vertical force, and 
its sideward horizontal ("medial") force. Figure || is based on 
the progressive force measurements. The transducer outputs 
were sampled and recorded at 2500 samples per second and 
were then averaged over 0.01 second intervals to produce, for 
each measured half-gait cycle, a single data set. Figure || is 
based on Bontrager's data set 3506a5 (data set a5 for person 
3506) ©j-. 

For Fig. g we modified the data set by adding, at the begin- 
ning, the last 9 points from the left-foot measurement; and at 
the end, the first 1 1 points from the right-foot measurement 
(under the plausible assumption that the unmeasured end of 
the previous left-foot gait cycle is the same as the measured 
cycle, and the unmeasured beginning of the next right-foot 
cycle is the same as the measured cycle). We have divided the 
full gait cycle into two half cycles A and B, as shown in Fig. 

Equation (Oh implies that the noise spectrum Sh{f) is pro- 
portional to \F(f)\/ / 4 , where F is the Fourier transform of 
F(t), the sum of the forces from the two feet (bottom panel of 
Fig. ||). As was discussed at the beginning of Sec. [TJ, the only 
features of F(t) that can contribute significantly at frequen- 
cies / ~ 10Hz are those that change on timescales r < 50 
msec. The net force F(t) varies only modestly on such short 
timescales, but its first time derivative dF/dt = F varies 
strongly: Shortly after placing his heel on the floor ("heel 
down"), the person begins to transfer weight from his trail- 
ing foot to his leading foot; this beginning of weight transfer 
entails a change AF ~ 4000N/s in the slope of the force 
curve (change of jerk) on a very short timescale r ~ 20 msec. 
A time St ~ 90 msec later, just before the trailing toe lifts off 
the ground ("toe off"), the weight transfer ends with a second 
sharp change of jerk — AF. 

These two jerk changes during each half gait cycle give 
rise to the following (approximate) form of the half-cycle's 
Fourier transform Fx/2 in the frequency band of interest, 
2.5Hz < / < 25Hz ~ 0.5/r: 



\ F l/2 



T 2sm(jr f6t)AF . 



(14) 
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TABLE I. Sudden changes of jerk, and ratio of vertical to hori- 
zontal force spectra for a half gait cycle, as inferred from Bontrager's 
data sets [01. 
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A numerical Fourier transform (FFT) of F(t) (Fig. reveals 
that other features with timescales > 5t = 90 msec pro- 
duce modulations of F(/) analogous to but no larger than 
the sin(7r/<5i). We shall ignore these modulations and cor- 
respondingly shall approximate s\n(nf8t) by its rms value, 
1 /V2, and therefore shall rewrite the above formula as 



/2l 



V2AF 



(15) 



Equation (|5j) is a fairly accurate representation of the half- 
gait-cycle Fourier transform, not just for the data set in Fig. 
pi but for all force-plate data that we have examined — data 
in the Biokinesiology literature [ |Ix||Tl"1 |, and unpublished data 
on two other subjects in Bontrager's laboratory JDj], 

In all of Bontrager's data sets except 3506a5 (Fig. [|), only 
one force plate was used rather than two, so the data sets show 
only the horizontal force produced on the floor by one foot 
during a half gait cycle, along with the times of heel down and 
toe off for the second foot as measured by a switch attached to 
the foot. To compute the force of the second foot, we assumed 
(in accord with Bontrager's advice) that its force history was 
the same as that measured for the first, but displaced in time 
as shown on the foot-switch recordings. We thereby modified 
each data set to include the force of the second foot, and com- 
puted the total force of the two feet (analog of segment A and 
segment B of Fig. [|). 

We have fit Eq. ( |l5| ) to the Fourier transforms of each of 
Bontrager's data sets for the total force of both feet during a 
half gait cycle. For each single-force-plate data set, and for 
the two half-gait-cycles (A and B) of dual-plate set 3506a5, 
we did a least-squares solution for AF = |# 1/2 1 (2tt/) 2 / ^2, 
and its rms fluctuations in the frequency range 2.5 to 20Hz. 
The results are shown in Table Q. 

Note that in each measured spectrum there are fluctuations 
of typical magnitude 30-60 per cent around the 1/ f 2 law of 
Eq. (|i"5|); and the fluctuations in the inferred AF from one gait 
cycle to another and from person to person are of order 30 per 
cent. If a person were to run rather than walk, the resulting 
AF might be larger, but presumably not by more than a factor 
~ 2; and running in the vicinity of a LIGO test mass should 



be much less common than walking. These variations in AF 
are modest contributors to our overall factor ~ 3 uncertainty 
in the gravity gradient noise. Based on Table [j| we shall use 
the value 



AF = 5500N/S 



(16) 



in our noise evaluations. 

For Bontrager's data sets we have also evaluated the ratio 
|-Fv|/|Fh| of the Fourier transforms of vertical force and hor- 
izontal force in the vicinity of 10 Hz; see the last column of 
Table |. The vertical spectra vary by a large factor from one 
person to another: male 3700 walks rather smoothly; female 
3506 strikes the floor sharply with her heel at heel down. Cor- 
respondingly, female 3506 produces a large change of force 
AF V at heel down on a short enough timescale to dominate the 
vertical spectrum, F v ~ (2tt f)^ 1 AF V ; whereas male 3700 
(and also male 3772) has a vertical spectrum dominated by a 
change of jerk and therefore falling off more rapidly with fre- 
quency, |F V | ~ (2tt/)" 2 AF v . As a result, female 3506 has 
a significantly larger ratio of vertical to horizontal |F| than 
the males. However, even for her large heel-induced |F V |, the 
resulting vertical contribution to the gravity gradient noise is 
smaller than the horizontal contribution — smalle r by the fac- 
tor - 0.1|F v |/|F h | - 0.5 discussed inSec.|lIAl. 



3. Noise spectrum 

For people walking in the vicinity of LIGO test masses, the 
sharp changes of jerk are not likely to occur in a periodic fash- 
ion to within a period accuracy of 0.05 sec, and correspond- 
ingly, in the vicinity of 10 Hz the jerks are not likely to su- 
perpose coherently. Therefore, we can approximate the sharp 
changes of jerk as constituting a random shot noise, for which 
the spectral density of the gravitational-wave noise from a sin- 
gle person will be 



% = 2 



F 



gait 



1/2 



\hx /2 \ = 



2 % /2aGAF 



(17) 



Here F ga j t is the gait-cycle period (about 1 second). In the 
first expression 2/F ga i t is the rate of (dual-jerk) "shots" (half 
gait cycles) for each of which hi/2 is the Fourier transform of 
h(t), and the second expression follows from Eqs. ( |l3| ) and 

For a number of walking people, each at a different distance 
r,; from the test mass and with a different angular location and 
direction of motion and corresponding factor a,, the noises 
will add in quadrature, producing 



Sh 



2y/2GAF 



1/2 



(18) 



Inserting the numerical values (discussed above) \a>i\ — 

M representative = V%, AF — 5500N/ S , L = 4km, F gait = 

Is, we obtain the noise spectrum ([[]) discussed in the introduc- 
tion and the abstract. 
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B. Motion of Floor and Ground 

Sharp changes in the horizontal force F(t) will produce 
deformations of the floor and ground, which become seismic 
waves as they propagate out through the earth. These defor- 
mations will produce gravity gradient noise correlated with 
that from the person's CofM. In this section we shall estimate 
this noise. 

We begin with the standard expression [0] for the gravi- 
tational potential at test-mass location xa, produced by the 
displacement ( of the floor and ground: 



* = -/ f^-ifw 

\X - X A \ 



G pC-dA 
ov \x — xa\ 



(19) 



Here —V • {pQ is the Eulerian change in density induced by 
the displacement (, the first integral is over the interior V of 
the floor and ground, and the second integral is over the sur- 
face layer of mass produced by ( on the surface dV of the 
floor and the adjoining ground outside the LIGO buildings. 

Integrating the first term by parts and cancelling the result- 
ing surface term against the second term of ( 19), we obtain 



G P Q 



1 



\x - xa 



d 3 j 



(20) 



The gravitational acceleration on test-mass A is minus the gra- 
dient of this with respect to xa- Since the only dependence of 
<£> on xa is through the combination x — xa, the gradient can 
be replaced by a derivative under the integral with respect to 
—x, thereby giving 



GpQ 



1 



xa\ 



(21) 



This is the dipolar floor/ground analog of Eq. (^|) for the 
gravitational acceleration produced by the person's CofM. 
By the procedure that led from (p|) to (|l0|), we ob- 
tain the gravitational-wave noise d 2 hfdt 2 produced by the 
floor/ground motion, to which we add the person's CofM 
noise (|l0|). Differentiating the result once in time, we obtain 



~dW 



G 



E 



mAi 



\x — xa\ 



d^j 



(22) 



Here (as should be obvious) in the first term x is the person's 
CofM location and in the second it is a location in the floor or 
ground. 

This equation expresses the gravitational noise in terms 
of sharp changes in the person's CofM momentum and the 
momentum density of the ground. By momentum conser- 
vation, any change of the CofM momentum M£j must be 
accompanied by an equal and opposite change of the total 



floor/ground momentum J v pQd 3 x. If the suddenly deposited 
momentum remains close to the person [within a distance 
■C r = (person's distance to nearest test mass)] during a 
time r ~ 1/(2/) — 50msec, then the gravity-gradient noise 
from the floor/ground will nearly cancel that from the person's 
CofM. If the deposited momentum spreads out over a distance 
^> r in the time r, then it will produce a neglible gravitational 
force on the test mass, and negligible gravity gradient noise. 

The deposited momentum moves outward through the floor 
and ground with seismic-wave speeds; it resides in a spread- 
ing, widening shell whose sharp outer edge moves at the seis- 
mic P-wave speed cp and fuzzy inner edge at a little less than 
the seismic S-wave speed cs (see the Appendix). 

The floor on which the person walks is a slab of reinforced 
concrete 20cm thick. In each corner and end station this slab 
begins about 6m from the test mass and extends outward to 
about 18m from the test mass, and transversely about 12m 
in each direction. In the corner station the slab begins about 
10m from the nearest test mass and extends on outward an ad- 
ditional 10 to several 10's of meters, in a complicated shape. 
The concrete has cp ~ 3700m/s and cs — 1700m/s; so in 
r ~ 50msec, the outer edge of the spreading shell could move 
a horizontal distance of 180m and the inner edge, 85m if the 
slab were that large. To the extent, then, that the deposited mo- 
mentum gets trapped in the slab for r ~ 50msec, it spreads 
through the whole slab; and since most of the slab is some- 
what farther from the test mass than the nearest person (about 
10 meters) and in somewhat different directions (off to the 
sides), the slab's sudden momentum change will produce a 
considerably weaker gravitational force on the test mass than 
is produced by the person's sudden momentum change. 

The momentum that passes through the thin floor and into 
the ground below spreads through the ground at much lower 
speeds than that confined to the floor: cs — 270m/s; cp ~ 
520m/s in Hanford's dry soils and 1700m/s in Livingston's 
water- saturated soils; cf. Tables II and IV of Ref. [|l[|. Corre- 
spondingly, in r = 50msec time, the inner edge of the spread- 
ing momentum shell travels a distance ~ 13m; and the outer 
edge, ~ 25m at Hanford and ~ 80m at Livingston. These 
distances, being comparable to, and much larger than the per- 
son's ~ 10m separation from the test mass, will cause the mo- 
mentum suddenly deposited in the ground to produce a some- 
what smaller gravitational noise than that of the person; there 
is no possibility for a strong cancellation. 

On the other hand, in some cases the Green's function for 
the momentum spread (cf. Appendix, and Figs. 2-4 of Ref. 
Jl4|]) is strongly localized near the inner edge of the spreading 
shell. In such cases, with the momentum having spread only 
a distance ~ 15m compared to the person's ~ 10m distance 
from the test mass, the ground's gravitational force on the test 
mass might be as much as half that of the person, thereby re- 
ducing the net noise by a factor 2 relative to that of the person 
alone. This is a significant contributor to our estimated factor 
3 uncertainty in the net noise. 
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C. Motion of Limbs 

Turn, next, to the noise produced by the person's quadrupo- 
lar gravitational field 



Quad 



3 GT 
2 



jk i i 
~ x 3 x k 



(23) 



[Eq. (§)], Here Ijk is the quadrupole moment about the per- 
son's moving CofM, and x'j and r' = \x'\ are the vector 
and distance from the CofM to a test mass. Replacing x'j by 
Xj — £j, where Xj reaches from the (fixed) CofM location at 
the center of a gait cycle to the test mass, and £j (t) is the dis- 
placement of the CofM relative to that fixed point, we obtain 



aad 



3 GljktijTik 3 GXjk 



2n^fe ) , 
(24) 



where £H is the component of £ along h = x/r, the direction 
to the test mass, and £rr is its projection orthogonal to h. 

This quadrupolar gravitational field produces interferome- 
ter noise via jerky changes of Tjk (noise "intrinsic" to the per- 
son's quadrupole moment) and via jerky changes of £ ("ex- 
trinsic" noise). 

The extrinsic noise is readily seen to be very small com- 
pared to that from the person's CofM gravitational field 

$CofM = -GM^/r 2 : 



cxt 
h 



< 9 Tjk 



gCofM 



2 Mr 2 



0.01 



/10m 



(25) 



Here we have used an obvious estimate of the person's 
quadrupole moment. 

The intrinsic noise, arising from jerky changes of Tjk, has 
contributions from both terms in Eq. (H). That from the sec- 
ond term is smaller by ~ 2£/r < 0. 15(T0m/ r) than that from 
the first term, so we shall ignore it. Taking the gradient of 
the first term to obtain the gravitational acceleration on the 
test mass, and proceeding as in the derivation of Eq. (|To|), we 
obtain the following expression for the quadrupolar noise in 
the interferometer as a sum over contributions from the test 
masses A — 1, 2, 3, 4: 



d 2 h ^ 3G / 



-lijUAiriAjTUAknAk 



(26) 



We choose the x axis along the progressive direction (direc- 
tion £ of motion), the y axis along the medial (transverse hor- 
izontal) direction, and the z axis vertically upward. Then, be- 
cause rriA and fiA are horizontal vectors, and the body's jerky 
motion is in the x-z plane, the noise ( p6[ ) arises solely from 
two components of the quadrupole moment, 

T-xx ~^j\^Ixx Izz) ) ^yy 7^(Ixx "F ^zz) ■ (27) 



Here Ijk is the second moment of the body's mass distribution 
(the integral of px'jx' k over the body). 

As the person walks, the dominant contributors to jerky 
changes of the quadrupole moment are the motions of his legs. 
(His arms are less massive and jerk less.) We divide each 
leg into two parts, the thigh (reaching from hip to knee) and 
the shank (reaching from knee to ankle); and we approximate 
each of these as a point mass located at its center of mass, 
thereby making an acceptably small error. Measurements dis- 
cussed below show that the quadrupolar noise at frequencies 
/ ~ 10Hz arises primarily from sudden (t ~ 0.5// ~ 50 
msec) changes of the thigh and shank accelerations Aafc each 
time the person's heel strikes the floor (heel down). The cor- 
responding sudden change of Ijk (second time derivative of 
Ijk) is 



AI 



jk 



2mV-Aaf 



2m t XjAa{ 



(28) 



where the superscripts s and t denote shank and thigh, m is the 
mass of shank or thigh, and x'j is the vector from the person's 
CofM to the center of mass of shank or thigh at heel down. 

Biokinesiologists measure the motions of thigh and shank 
in two ways: by videotaping markers placed on them (position 
measurements), and via accelerometers placed on them (ac- 
celeration measurements); for a pedagogical discussion see 
[|T3|]. Because of noise introduced when taking time deriva- 
tives of the position data, the position measurements cannot 
give reliable measures of acceleration on the short timescales 
r < 50ms of concern to us Jl^]. Therefore, for the thigh and 
shank accelerations we have relied on accelerometer measure- 
ments as reported by Ge Wu in Fig. 16-11 of Ref. JT3|]. It is 
obvious from that figure that the dominant contributions to 
the Fourier transform a b k of a h k (t) (for k = x,y,z, b = t, s) 
arise from sharp changes at heel down. We have Fourier trans- 
formed a\ (t) and found that, to within a factor ~ 2 over the 
range 2.5Hz < / < 25Hz, \a b k \ oc 1//. This implies that, to 
adequate accuracy for our purposes (factor 2) we can regard 
the quadrupolar noise as due to sudden changes of accelera- 
tion I Aotl = 2Trf\a b k \ at heel down. 

Table || shows the values of \Aa k \ at heel down for Wu's 
typical 60kg individual, as inferred from our Fourier trans- 
forms of her Fig. 16-11; it also shows the values x'\ of the 
center of mass location of leg and shank at heel down for a 
typical individual, as given in Appendix A of Ref. [Q.Q We 
expect these numbers to vary, from one adult individual to an- 
other, by no more than a factor ~ 2. 

By inserting the numbers from Table |li| into Eq. ( p8| ) and 
thence into Eq. (|27|), and adding terms in absolute value so 
as to get upper bounds, we obtain for the sudden (timescale 
< 50 ms) heel-down changes of the second time derivative of 
the person's quadrupole moment: 



*Beware: Biokinesiologists (influenced by the Biomechanics liter- 
ature) use different axis conventions from physicists: y and z are 
interchanged so their y is vertical and z is medial. 
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TABLE II. Thigh and shank properties at heel down. Masses m 
(in kg) and positions x' , z' relative to the body's center of mass (in m) 
were taken from Appendix A of pip. Changes \Aa,j | of acceleration 
(in m/s 2 ), on the timescales r 5- 50 ms, were computed from Fourier 
transforms of Fig. 1 1-16 of [113]. 





m 


x' 


z' 


\Aa x \ \Aa z \ 


Thigh 


5.7 


0.08 


-0.3 


15 6 


Shank 


3.5 


0.2 


-0.6 


11 6 


\ax xx 


< 35kg 


mV 2 , 


\Ai yv \ 


< 25kg m 2 s- 2 . (29) 



By then Fourier transforming Eq. (pq), we obtain the fol- 
lowing upper limit on the quadrupolar noise | h Aa | induced by 
the heel-down changes of thigh/shank acceleration: 



< 



9G\A1 XX \ 
2(2vr/) 5 J Lr 4 



(30) 



where r is the distance to the nearest test mass. 

During half a gait cycle (one heel down), the Fourier trans- 
form of the person's CofM-induced noise is [Eqs. ( |l3| ) and 



\h G 



oiM 



2G|AF| 
Lr 3 {2irfY 



(31) 



where we have set the angle-dependent factor \a\ to its repre- 
sentative value, V2- Near 10 Hz, the ratio of the quadrupolar 
noise (BOb to this CofM noise is 



h Aa \ ^9 27rf\AX xx 



\hco 



I'M 



< _ 

4 



\AF\ 



1 

10 



10m \ / / 



V 10Hz 



(32) 



This is significantly less than one independent of the factor ~ 
2 errors and variabilities of both noises. Thus, as was asserted 
in the Introduction, the dominant noise is caused by jerkiness 
in the person's CofM motion. 



III. DOORS, FISTS, AND VEHICLES 

We turn, now, to gravity gradient noise produced by the im- 
pulsive stopping of a massive, horizontally moving object — 
most especially a slamming door, a fist striking a wall, or a 
stopping automobile. (The impulsive stopping of vertical mo- 
tion produces a much weaker signal than horizontal motion; 
cf. the discussion following Eq. (10).) Our analysis is a vari- 
ant of that originally given by Spero [g] and reaches the same 
conclusions. 

Since these impulsive events are not likely to occur repet- 
itively and continually (by contrast with the gait cycles of 
human walking), they are more appropriately analyzed as a 
single impulsive gravitational-wave signal than as a stochas- 
tic noise. The standard formula for the amplitude signal to 
noise ratio S/N produced by such an impulsive signal h(t) in 



a LIGO interferometer is (e.g., Eq. (29) of Jl 5p with a factor 
2 correction) 



N 2 



A\h\ 2 



df, 



(33) 



where Sh is the one-sided spectral density of the interferom- 
eter's total noise and h is the Fourier transform of the signal 
h(t). 

The signal is given by 



(34) 



[Eq. (|TT|)], where a is the same angle-dependent factor as we 
met for a person's CofM motion [Eq. (|l2|)], M is the object's 
mass, £ is its displacement while stopping, and r is its distance 
mass from the nearest interferometer test mass. 

For a slamming door or a fist striking a wall, it is the ve- 
locity v = £ that changes suddenly, by some amount Av. 
Correspondingly, the Fourier transform of the signal h(t) is 



\h\ 



GM\aAv\ 
J Lr 3 (2 7 r/) 4 



(35) 



For the benchmark "advanced" LIGO interferometer, 
we can approximate the noise curve (Fig. [j]) by Sh = 
So[(fo/ff + (fo/f) 20 }, where S - 10- 45 /Hz and f = 
10Hz. Inserting this Sh and expression (|5|) into Eq. ([33|) and 
integrating, we obtain 



1.2GM|aAu| 



2^LrH2^f f^Sj2 
MAv \ /10m\ 3 
5kg m/s J \ r J 



(36) 



Here we have used our representative value \[2 for \a\. This 
is the noise level discussed in the introduction. 

Turn to automobiles and other vehicles. Under normal 
(non-collisional) motion, the velocity of a vehicle cannot 
change significantly on a timescale of 50 ms. The accelera- 
tion, however, can so change, and will change most strongly 
when the vehicle comes to a stop, e.g. when parking. 

The sudden change Aa = A£ of acceleration, when the ve- 
hicle comes to rest, will produce the following Fourier trans- 
form of h(t): 



\h\ 



GM\aAa\ 
Lr 3 (27r/) 5 



(37) 



Inserting this \h\ and the advanced-interferometer Sh(f) into 
Eq. (33 ) and integrating, we obtain 



!3_ r 



GM\aAa\ 



2^Lr3(2^/ ) 4 V^7o 
( M \ (\Aa\\ /30m 
0£g) \T~ 



if I 

2tonne / 



(38) 
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Here we have used our representative value \[2 for \a\, and g 
is the acceleration of gravity. This is the noise level discussed 
in the Introduction. 

The gravitational signal from a slamming door, striking fist, 
or stopping vehicle will be mitigated to some modest extent 
by an opposite signal produced by the momentum deposited 
in the "reaction mass" (the wall, floor, and/or ground). How- 
ever, as for human walking, the deposited momentum spreads 
over such a large spatial region in a time 1/ 2/ ~ 50ms, that 



the mitigation will not be significant; cf. Sec. [I B and the Ap- 
pendix. 



IV. CONCLUSIONS 

In this paper we have identified what we believe to be the 
dominant gravity gradient noise due to normal human activi- 
ties; we have estimated its spectrum and its strength to within 
accuracies of a factor ~ 3; and we have discussed the impli- 
cations of this noise for the size of the human exclusion zones 
around LIGO's test masses in the era, ca. 2010, of "advanced" 
interferometers. Until that era, human gravity gradient noise 
is not likely to be a serious issue for LIGO. 

Our formulas and estimates can provide a basis for the de- 
sign of the facilities of other earth-based gravitational-wave 
detectors. 
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APPENDIX A: GREEN'S FUNCTION FOR SPREADING 
MOMENTUM IN FLOOR AND GROUND 

When a time varying force Fj (t) is applied to the surface of 
the earth at a location of, the time varying displacement that it 
produces in the ground at a location x is given by 



Ci(x, t) = Q- t j 9" & t - t'; f )F j {t')dt l 



(Al) 



Here gfj is the elastodynamic Green's function for a unit-step- 
function (Heaviside-function) force, F(t') = H(t')ej. Seis- 
mologists focus on gfj rather than on the physicists' usual 



delta-function-sourced Green's function gij = dgfj /dt be- 
cause gfjS Heaviside steps at the various seismic propaga- 
tion fronts are more easily visualized and compared with each 
other than g^-'s delta-function spikes. 

The Green's function g*j for a homogeneous medium ("ho- 
mogeneous half space") has been computed analytically by 
Lane Johnson [ff3], up to a complicated integral; and Johnson 
has evaluated it numerically for several representative geome- 
tries; see his Figs. 2-A. Chao jl6| ] has derived expressions for 
gfj when both force point and field point, x' and x, are at the 
surface of a homogeneous half space; and Ma and Huang jl7| ] 
have computed gfj for layered media. 

Regardless of the nature of the medium, momentum con- 
servation requires that 



Pi® 



d 
dt 



n^d 3 x 
9 dt 

d?_ 

dt 3 



P9 



Z&t-t'^^x 



F k {t')dt' , 
(A2) 



where V is the entire volume of the medium. Since this must 
be true for every applied force, it must be that 



d?_ 

dt 3 



pgg(x, t - t'; x')d 3 x = S(t - t')6 jk 



(A3) 



Causality requires that gfj vanish everywhere for t < t' and 
be nonzero for x arbitrarily near x' when t — t' is arbitrarily 
small but positive; these facts, combined with Eq. ( A3 ) imply 



it - t') 2 

pgfj(x, t - t'; x)d 3 x = ± — ^-Hit - • (A4) 



Now consider the noise h(t) produced in a gravitational- 
wave interferometer by a walking person, whose feet at 
location x' produce a horizontal force Fj(t) on the floor 
and thence on the ground beneath the floor. By (i) taking 
three time derivatives of Eq. ([22]), with the foor and ground 
displacement £ expressed as an integral over the elastody- 
namic Green's function [Eq. (|A1[)], (ii) using force balance 
d(M£j)/dt = —Fj for the floor and person, and (iii) set- 
ting d 2 Fi{t)/dt 2 — AFiS(t), where t = is a time of sharp 
change of jerk at the beg inning or end of the walking person's 
weight transfer (cf. Sec. Q A 2), we obtain the following: 



d 6 h . ■ J ,, i . ( I 

7 rriAiAFk < o{t)d jk 



dt (i 



dt 3 



pgf k (x,t;x") 



\x - xa\ 



X - XA\ 



d 3 3 



(A5) 



This equation exhibits the momentum-flow features dis- 
cussed in the text following Eq. (j22|): The third time derivative 
of the Green's function g^ k (x,t\x') is significantly nonzero 
only in the expanding shell discussed in the text — a shell 
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whose sharp outer edge travels at speed cp and fuzzy inner 
edge a bit slower than cs; cf. Figs 2-4 of Johnson Jl4|]. Cor- 
respondingly, the contribution of the floor and earth to h is 
confined to that shell. When that shell is small compared to 
the separation \x' ~ xa\ between the person and the nearest 
test mass A, the double gradients in Eq. (A5) are nearly equal, 
and momentum conservation as embodied in Eq. ( A3 ) guar- 
antees that the two terms in (A.5) (the person noise and the 
floor & ground noise) will nearly cancel. When the shell is 
comparable in size to the separation \x' — xa\, the two terms 
will cancel partially but not strongly. When the shell is large 
compared to \x' — Xa\, the second term (floor & ground noise) 
will be negligible compared to the first (person noise). 

Equation (|A5|), together with the explicit expressions for the 
Green's functions in Refs. [ |l4|jr^ , [l7| ], could be used to com- 
pute quantitatively the partial cancellation of person noise and 
floor & ground noise. We have not done so, since uncertain- 
ties elsewhere in our modeling are comparable to or larger 
than the errors in the above rough estimates. 
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